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(54) Methode pour gene re r un maillage 3D 

(57) -La methode peut s'appliquer par exemple a 
une structure geologique definie a partir de don- 
nees disponibles acquises soit par des mesures, 
soit par des analyses et des interpretations de don- 
nees sismiques ou des diagraphies realisees dans 
despuits. Ellecomporte une moderation prealable 
des limites (bords de la zone, discontinues ou 
failles par exemple) par des surfaces triangulees, 
enveloppant au mieux ces limites, et un fractionne- 
ment de la zone d'etude en blocs ayant au plus six 
faces, dont les faces correspondent a des disconti- 
nuites ou des failles principales de la zone preala- 
blement triangulees. Chaque bloc est ensuite sub- 
divise en un ensemble de cellules elementaires ne- 
cessaires pour la modelisation de la zone par ele- 
ments ou differences finis, en realisant un maillage 



de chacune des faces du bloc par interpolation a 
partir de ses bords de facon a obtenir une grille, sui- 
vi d'un peuplement en noeuds de chaque bloc dans 
son epaisseur, par interpolation a partir des noeuds 
des differentes grilles formees sur les differentes fa- 
ces, ces operations etant accompagnees de proce- 
dures particulieres d'homogeneisation pour reajus- 
ter la taille des cellules a I'interieur de chaque bloc 
suivant des regies definies. Les blocs mailles sont 
ensuite reassembles. 

Application notamment a I'etude de gisements pe- 
troliferes. 
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Description 

La pr6sente invention concerne une methode pour generer un maillage 3D dans le but de modeliser un corps 
quelconque dont ia geometrie est connue et imposee a priori tel que par exemple une structure geologique en tenant 
5 compte fidelement de sa geometrie dans toute sa complexite naturelle. Le maillage realise par la methode s'accorde 
bien dans le cas d'une structure geologique, avec la position des interfaces entre les couches, avec les failles meme 
de forme complexe, les biseaux et autres irregularites de structure, en les representant avec leur conformation reelle. 

La methode selon I'invention trouve des applications notamment pour la modelisation en geomecanique, pour la 
caracterisation des reservoirs : couches de sediments, stratigraphie, et aussi pour la simulation des ecoulements dans 
10 les bassins et reservoirs. 

I) ETAT DE LA TECHNIQUE 

La modelisation d'une structure geologique est realisee a partir d'un ensemble de donnees disponibles acquises 

15 soit par des mesures, soit par des analyses et des interpretations de donnees sismiques ou des diagraphies realisees 
dans des puits. Ces donnees se presentent sous la forme d'un nuage de points entre lesquels, au mieux, il convient 
d'installer un maillage. La formation d'un maillage en 3D joue un role crucial pour faire des simulations numeriques 
dans differentes applications geologiques. Beaucoup de recherches ont ete menses sur le sujet principalement en 
ingenierie de reservoir en raison de ses applications dans le domaine petrolier. 

20 On connait de nombreuses techniques de maillage 3D pour des simulations numeriques par elements ou diffe- 

rences finis et des ensembles de logiciels de simulation crees pour leur mise en oeuvre qui sont adaptes a des appli- 
cations et respectent des contraintes determinees. 

II est connu pour certaines applications d'utiliser des maillages de type tetraedrique non structure. II est connu 
aussi de faire du maillage hexaedrique de peau ou de surface et de repartir des points dans le maillage soit par 

25 interpolation transfinie, soit encore par resolution d'equations aux derivees partielles. Les surfaces sont g6n6ralement 
creees a partir de fonctions parametriques et les points du maillage sont obtenus simplement par echantillonnage de 
ces fonctions. Une telle methode ne convient pas gen£ralement pour des applications dans les g£osciences ou les 
volumes a modeliser ont des formes le plus souvent tres complexes, difficiles a traduire par des fonctions parametriques 
compte-tenu de la precision qui est recherchee. Aussi prefere ton approcher les formes complexes par des elements 

30 de surface triangules. 

Differentes techniques de maillage appliquees a la modelisation de reservoir, sont decrites par exemple dans : 

Georges RL : "Generation Automatique de Maillages - , Editions Masson, 1991; 

35 - Guerillot D.R. et al: An ineractive 3D Mesh Builder for Fluid Flow Reservoir Simulation; in SPE - Computer Appli- 
cations, Dec. 1 993; 

Aziz K.; Reservoir Simulation Grids : Opportunities and problems. SPE-JPT, July 1993; ou 

40 - Heineman Z.E. et al : Gridding Techniques in Reservoir Simulation. Intl Forum on Reservoir Simulation. Alpbach 
Austria, Sept. 12-16, 1988. 

Cependant, la plupart des techniques de maillage connues, et de la les simulations qui en resultent, ne conviennent 
bien que pour modeliser des structures geologiques relativement simples. Les irregularites geologiques courantes 

45 dans les zones-reservoir modeiisees, telles que les couches ayant de grandes variations d'epaisseur, les failles ou les 
biseaux, sont rarement pris en compte ou alors seulement sous une forme grandement simplifiee. Dans les modeles 
mailies obtenus dans le cadre des geosciences, les failles apparaissent verticalement; au mieux, elles sont inclinees 
et planes. Or il est important de rendre fidelement, la geometrie en 3D des sediments car ils possedent des caracte- 
ristiques physiques propres, et aussi la geometrie des failles et autres discontinuites mecaniques. 

50 La methode selon I'invention permet de generer un maillage 3D respectant la geometrie d'une zone d'un milieu 

complexe dont la geometrie est connue et imposee a priori, tel qu'une structure geologique, dans le but de realiser un 
modele representatif de ce milieu. La structure est .definie a partir de donnees disponibles acquises soit par des me- 
sures, soit par des analyses et des interpretations de donn6es sismiques ou des diagraphies r6alis6es dans des puits. 
Elle comporte une modelisation prealable des limites (bords de la zone, discontinuites ou failles par exemple) par des 

55 surfaces trianguiees, enveloppant au mieux ces limites, et un f ractionnement de la zone d'etude suivant en blocs ayant 
au plus six faces, dont les faces correspondent a des discontinuites ou des failles principals de la zone pr6alablement 
trianguiees. 

Elle est caracterisee en ce que'elle comporte : 
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la subdivision de chaque bloc en un ensemble de cellules elemental res necessaires pour modeiiser la zone par 
elements finis, par une reticulation en 3D de chacun des difterents blocs obtenus, subdivision qui comprend : 

a) le maillage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une grille 
5 intermediate, suivi d'une projection sur la surface trianguiee des noeuds de cette grille intermediate suivant 

une direction definie (perpendiculairement a elle-meme par exemple); 

b) le peuplement en noeuds de grille de chaque bloc dans son 6paisseur, par une methode d'interpolation a 
partir des noeuds des differentes grilles formees sur les differentes faces triangulees de chaque bloc, et 

10 

c) un rearrangement dans le positionnement des noeuds des cellules eiementaires crepes, obtenu (de facon 
iterative, generalement) par relaxation dans le positionnement des noeuds, en fonction de criteres predefinis 
tels que la r£gularite, le choix de la taille relative des differentes cellules et le mode de progression de leurs 
tailles, la direction de projection etc. 

15 

Chaque bloc etant ainsi reticule en 3D, on peut proceder a un reassemblage ou reempilement de tous les blocs. 

On peut proceder successivement a une premiere operation de relaxation dans le positionnement des noeuds sur 
chacune des differentes faces triangulees, le peuplement en noeuds internes etant effectue a partir des grilles relaxees, 
et ensuite a une deuxieme operation de relaxation dans le positionnement de ces noeuds internes. 
20 De preference, la premiere operation de relaxation est effectuee par des phases iteratives de relaxation-projection. 

La methode selon ['invention permet d'obtenir un modele qui respecte tres precisement la geometrie souvent 
complexe de structures geoiogiques telles que des champs petrolifdres, ce qui a pour effet que les simulations qui 
peuvent etre faites a partir de ce modele rendent plus fidelement compte des phenomenes. La methode utilis6e permet 
en outre une economie tr£s importante de temps de calcul pour realiser les operations de maillage 3D habituellement 
25 tres longues avec les methodes existantes. 

D'autres caracteristiques et avantages de la methode selon Pinvention, apparaitront a la lecture de la description 
ci-apres de modes de realisation d6crits a titre d'exemples non limitatifs, en se ref6rant aux dessins annexes ou : 

les Fig. 1 et 2 montrent des exemples de relaxation-projection sur un maillage surfacique 4X4; 

30 

la Fig. 3 montre un cube deforme avec des faces triangulees; 

la Fig. 4 montre I'enveloppe externe d'un cube maill6; 

35 - les Fig.5 et 6 montrent les surfaces internes d'un cube mailie respectivement avant et apres relaxation; 

la Fig. 7 montre des courbes relevees sur cinq coupes d'une structure geologique de tres grandes dimensions, sur 
laquelle la methode selon I'invention a ete testee, qui constituent une partie des donnees initiales dont on dispose, 
sur cette structure; 

40 

la Fig.8 montre d'autres donnees initiales sur la mdme structure, en I'occurrence, des cartes d'isobathes et les 
emplacements de 4 puits; 

la Fig.9 montre une moderation surfacique ( a base de surfaces triangulees) de la structure geologique etudiee; 

45 

la Fig. 10 montre un exemple de decoupage de la structure en blocs; 

la Fig. 11 montre le resultat d'un maillage bloc par bloc de I'enveloppe externe de la structure geologique; 

so - la Fig. 12 montre le resultat d'un maillage bloc par bloc des surfaces internes de la structure geologique; 

la Fig.1 3 illustre a titre d'exemple comment on prend en compte un biseau dans le maillage de la structure geo- 
logique; 

55 - la Fig. 1 4 illustre a titre d'exemple comment on prend en compte une faille dans le maillage de la structure geolo- 
gique; 

la Fig. 15 illustre la methode d'interpolation 2D de Coons; 
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la Fig. 16 montre les Elements pris en compte dans la methode d' interpolation 3D de Coons; et 

la Fig. 17 illustre le mode de relaxation-projection des noeuds de maille. 

s La structure geologique de grandes dimensions illustree par les Fig.7-10, correspond a un dome anticlinal com- 

portant six interfaces entre des sediments. 

La realisation d'un maillage en 3D par elements finis d'une zone a partir d'un ensemble de donnees initiales con- 
nues pour la zone : interfaces entre differentes couches ou horizons, failles principales, carte en courbes de niveau 
avec les emplacements de puits passe par une succession d'etapes qui seront decrites ci-apres. 

10 

I) Moderation des limites ou bords de la zone etudiee et de ses principales discontinues. 

On realise une modelisation des surfaces : bords, interfaces, failles en utilisant une logictel specialise connu tel 
que GOCAD (marque d6posee) sous la forme d'une juxtaposition continue d'6iements de surface triangulares. Dans 

'5 un premier temps on forme une ebauche qui s'adapte au mieux aux donnees initiales. De maniere iterative, on deplace 
les sommets de la surface en utilisant un logiciel d'interpolation et on densifie localement les Elements de surface 
triangulaires, de facon a plaquer au mieux la surface sur le semis de points (donnees d'entree). Les elements de 
surface ainsi cr66s ne se raccordent generalement pas comme ils le devraient : les interfaces ne sont pas limitees aux 
intersections et au bord du domaine d'etude et elles peuvent se croiser en debordant. On calcule alors les lignes 

20 d'interception entre surfaces et on limite celles-ci aux courbes d'intersection. La Fig 7 montre le r6sultat final de cette 
modelisation des limites de la structure et Ton peut remarquer que la plupart des triangles des surfaces triangulees du 
biseau sont localises dans les regions de forte courbure. 

A ce stade, il n'est pas encore possible de produire un maillage 3D tenant compte des contraintes geologiques et 
numeriques. Pour y parvenir, on precede a un fractionnement de la zone. 

25 

II) Fractionnement de la zone etudiee 

On decoupe d'abord la zone en macro-blocs ayant au plus six faces. II s'agit generalement de blocs de forme 
hexaedrique (a six faces), mais le nombre de faces peut etre inferieur a six et descendre jusqu'a deux. Dans ce cas, 

30 pour les besoins de la methode, le bloc est quand mdme considered comme un block hexaedrique mais ayant des 
faces degenerees (reduites a un point ou a une courbe). Le mode de decoupage est generalement suggere par la 
geometrie des horizons, des principales failles et leur interaction. Pour constituer ces macro-blocs 3D, on a besoin 
parfois d'introduire des f rontidres artificielles. C'est n6cessaire par exemple quand un bloc a plus de six faces naturelles 
et ne satisfait done pas a la regie choisie d'etre au plus hexahedrique. Dans un tel cas, on reunit deux faces en une 

35 face complexe pour satisfaire le critere, ou bien on introduit des surfaces de coupure artificielles. 

Un exemple de decoupage de la structure geologique etudiee en macro-blocs est montre a la Fig. 10, ou Ton a 
extrapole la gaille jusqu'a la frontiere de la structure et verticalement projete I'intersection du biseau sur les autres 
surfaces. 

40 III) Maillage 

On realise une reticulation ou maillage en 3D separement de chacun de ces macro-blocs de facon a le subdiviser 
en blocs eiementaires qui vont former les cellules pour la modelisation par elements finis recherchees. Cette operation 
est realisee en deux temps. 

45 

1) Maillage des faces 

a) on cree tout d'abord une grille sur chacune des cinq ou six faces de chaque macro-bloc. Pour cela, on forme 
une grille intermediaire par interpolation a partir des quatre bords de la face, en utilisant la methode d'interpolation 
50 transfinie de Coons dont le principe est rappeie ci-apr£s. 

Etant donne quatre courbes parametriques f1(u), f2(u), g1(v), g2(v) (0 2 u, v 2 1 ) d6finissant quatre bords jointifs 
(Fig. 15), on introduit une formule pour calculer la surface la plus tendue qui passe par les quatre bords (qui interpole 
ces bords) : 

55 

S(u,v) = 0-u)g-\(v)+ug2{v)+V~v)n(u)+vf2(u)-[V-uW-v)P{0,0) 

+(1 -u) vP(0,1 )+u(1 -v) P(1 ,0)+uvP(1 ,1 )] (1 ) 
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Les sommets P(i, j) sont les quatre coins de la surface. 

On peut alors obtenir un maillage regulier NxM de cette surface par simple echantillonnage en prenant pour som- 
mets les points S(i, j) correspondant aux ui = i/(N-1 ), vj = j/(M-1 ), i=0 a N-1 , j=0 a M-1 . La formule 1 restreinte au points 
du maillage devient : 

5 

SU.J) = (1-^)510)^920^1-^)^1(0+^ /2(/) 

Dans le cas illustrS, les courbes f1, f2, g1, g2, sont obtenues a partir du logiciel GOCAD prdcite et sont done des 
lignes brisees. Pour obtenir un maillage regulier MxN, on Schantillonne r^gulierement (en abscisse curviligne ou lon- 

15 gueur de segment) les courbes f 1 et f2 en N points, et les courbes g1 , g2 en M points, puis on applique la formule (2). 
On peut tout aussi bien obtenir un maillage a progression geometrique sur chaque direction, simplement en imposant 
cette contrainte au niveau du reechantillonnage des courbes. 

Ce maillage intermSdiaire sauf cas exceptionnel ne convient pas car il n'y a pas de raison qu'il coincide a ce stade 
avec la surface qui a pr^alablement ete triangulee par le logiciel de surfacage. L'ajustement n6cessaire consiste a 

20 projeter les noeuds de la grille intermediate sur la surface triangulee orthogonalement a celle-ci, ou suivant d'autres 
directions , de facon a obtenir une nouvelle grille avec des noeuds places sur la surface triangulee. Quand les quatre 
courbes frontieres sont regulierement echantillonnees, la methode d'interpolation de Coons presente I'avantage de 
produire une grille reguliere en ce sens que sur chaque courbe isoparametrique tous les segments echantillonnes ont 
la meme longueur. Plus generalement la methode de Coons conserve les rapports de taille entre elements cons6cutifs 

25 dans une direction donnees. Ceci permet d'avoir par exemple des maillages r6guliers ou respectant certaines pro- 
gressions dans la taille des mailles. Cependant la grille obtenue par cette methode connue perd sa regularity (ou le 
critere recherchy) une fois qu'elle est projet^e sur la surface triangulee. 

2) Rearrangement des noeuds de face 

30 

Pour retrouver le critdre desire dans la distribution des noeuds de la grille projetee, une regularity satisfaisante 
par exemple, il faut deplacer les noeuds du maillage sur la surface triangulee. Pour cefaire, on procede a une relaxation- 
projection de la grille par ajustements successifs comme indique ci-apres (dans le cas ou le critere est la regularity), 
en relation avec les Fig.1 , 2. 

35 Etant donne un bloc maille, le but cette relaxation est de rendre son maillage le plus regulier possible au sens ou 

n'importe quelle courbe isoparametrique du maillage sort regulierement echantillonnee. Pour ce faire, on deplace de 
facon iterative tout point du maillage de maniere a le positionner au mieux a 6gale distance de ses deux voisins sur 
chacune des trois directions (u.v.w) de voisinage. Soit P le point qu'on cherche a deplacer. Pour chaque direction i, (i 
= u, v,w) (Fig. 17), on trouve le point Pi qui se situe a egale distance des deux voisins de P sur la courbe de direction 

40 j, (j = u, v,w). Le point P est alors remplac6 par le barycentre des trois points Pi (i = u, v, w). 

Apres avoir deplac6 tous les points du maillage de la sorte, ceux-ci ne sont plus sur la surface triangulee. On les 
projette alors a nouveau sur la surface triangulee et on repete le processus de relaxation -projection jusqu'a atteindre 
la regularity souhaitee. 

Ceci est illustry sur I'exemple des Fig. 1, 2, correspondant a une face du cube deform6 de la Fig. 3 qui est un 
45 exemple synthetique permettant de bien illustrer les diff£rentes etapes de la methode. Cette surface a 6te obtenue 
par simple deformation d'une face carree en preservant ses bords ( le bord de la surface reste un carry). Sur ces 
figures, on peut voir trois maillages 4X4. Le premier est un maillage plat issu de I'interpolation de Coons appliqu6e 
aux quatre bords (formant un carrd). Le deuxieme est le maillage fin dont les noeuds sont obtenus par projection de 
ceux du premier maillage. On constate que ce maillage a perdu de sa regularity. Le troisiyme est le maillage regulier 
50 tracy en gras et dont les noeuds sont situys sur la surface. II a ety obtenu a partir du deuxiyme maillage aprys une 
syrie de projections et relaxations. La Fig. 4 montre le rysultat de ces operations appliquyes a I'ensemble des faces 
du cube dyforme. 

3) Peuplement en noeuds de grille de chaque bloc dans son ypaisseur : 

55 

Pour prolonger en profondeur les grilles formyes sur les faces des blocs (Fig.5), on utilise une gynyralisation au 
volume de la mythode d'interpolation transfinie prycedente. 

Etant donny 12 courbes paramytriques f 1(u), f2(u ), f3(u), f4(u), g1 (v), g2(v), g3(v), g4(v), hi (w), h2(w), h3(w), h4 
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(w) (0 2 u, v, w 2 1) definissant les 12 bords d'un bloc tel que montre dans la Fig.16, une generalisation de la formule 
de Coons fournit une parametrisation volumique a I'interieur du bloc. La formule de Coons dans le cas volumique est 
la suivante : 

V( u,v, w) = ( 1 - v)(1 - w) fl ( u)+( 1 - v) wf2( u)+ vwf3( u)+ v(1-w) f4{ u) 
+{^-u)(^-w)g^{v)+(^-u)wg2(v)+uwg3(v)+u(^-w)g4(v) 
10 +^-u){1-v)M(w)+0'U)vh2(w)+uvh3(w)+u^-v)hA{v) 

+C(u,v,w) (3) 

avec 

75 

) C(u,v f *v)=-2[(1-u)(1-^ 

+(1-ty)v(1-w)P(0 > 1 l O)+(1-iy)vwp(O l 1,1)+t;(1-v)(1-w)P(1 I 0 > 0) 

20 

+u(1 - v) wP{\ ,0, 1 )+uv(1 - w) P(1 ,1 ,0)+uvwP(1 , 1 , 1 )] 
Les sommets P(i, j, k) ( i, j, k= 0 ou 1 sur la Fig 16) sont les 8 coins du bloc. 

En combinant les equations (1) et (3) on aboutit a une expression de la formule de Coons dans le cas volumique 
25 en fonction des 6 surfaces de bord du bloc, S1 (u, w), S2(u,w), S3(v,w), S4(v,w), S5(u,v), S6(u,v). Mais ici, les surfaces 
sont aussi suppos6es etre obtenues par interpolation de Coons a partir de leurs courbes de bord. La formule est la 
suivante : 

30 V(i/,v,mO = |[(1-WS1(u,mO+^ 

+( 1 - w) S5( u, v)+ wS6( u,v)]+ 1 C{ u, v, w) (4) 

Un echantillonnage quelconque de V(u,v,w) fournit alors un maillage EF du bloc. 
35 Cette interpolation en profondeur donnerait une distribution reguliere de noeuds si elle etait appliqu^e a partir des 

noeuds des grilles formes sur les differentes faces du bloc. Ce n'est pas le cas et la distribution obtenue est g6n6ra- 
lement irr6guliere car la mSthode d'interpolation en profondeur est appliquSe aux noeuds reports par relaxation sur 
les surfaces triangulees, comme decrit dans I'etape 2 precedents. 

40 4) Rearrangement des noeuds de bloc 

Aussi, pour retrouver une r6gularite suffisante (ou respecter un autre critere choisi), la methode selon I'invention 
comporte egalement une relaxation de la grille 3D de la maniere indiquee pr6c6demment. Les Fig.5, 6 montrent le 
maillage interne du cube deform^ avant et apres relaxation. Le maillage de la Fig.6, on le note, possede une regularity 
45 bien meilleure que le maillage de la Fig.5. 

Quand chacun des blocs precedents de la structure etudi6e a 6X6 convenablement mailie en volume, on reassem- 
ble tous ces macro-blocs apres leur reticulation 3D. 

La methode qui vient d'etre ddcrite peut etre utilisee pour mailler des h£terogeneites diverses telles qu'un chenal 
ou encore pour mailler radialement un puits deVie\ 
50 On a decrit des modes de realisation de ia methode appliquee au maillage d'une structure geologique. On ne 

sortirait pas du cadre de I'invention, cependant, en utilisant la methode pour le maillage de tout corps ou structure dont 
la g6om6trie est connue et impos6e a priori. 



55 Revendications 

1. Methode pour g6n6rer un maillage 3D dans une structure dont la geometrie est connue et impos6e a priori telle 
que par exemple une structure geologique definie a partir de donnees disponibles acquises par des mesures, par 
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des analyses et/ou des interpretations de donnees sismiques ou des diagraphies realisees dans des puits au 
travers de la structure, dans le but de realiser un modele representatif de cette structure, comportant une mode- 
lisation prealable des limites par des surfaces triangulees, enveloppant au mieux ces limites, et un fractionnement 
de la zone en blocs ayant au plus six faces, dont les faces correspondent a des discontinuites ou des failles 
principales de la zone prealablement triangulees, la methode etant caractensSe en ce qu'elle comporte : 

la subdivision de chaque bloc en un ensemble de cellules elementaires necessaires pour modeliser la zone 
par elements finis, par une reticulation ou maillage en 3D de chacun des differents blocs obtenus, subdivision 
qui comprend : 

a) le maillage de chacune des faces du bloc par interpolation a partir de ses bords, de facon a creer une 
grille intermediate, suivi d'une projection sur la surface triangulee des noeuds de cette grille intermediaire 
suivant une direction definie, de facon a generer des grilles de face; 

b) le peuplement de chaque bloc en noeuds internes par interpolation a partir des noeuds des differentes 
grilles formees sur les differentes faces triangulees de chaque bloc, de facon a obtenir une distribution 
reguliere des noeuds dans le volume de chaque bloc; et 

c) un rearrangement dans le positionnement des noeuds des cellules elementaires creees, obtenu par 
relaxation dans le positionnement des noeuds, en fonction de criteres predefinis. 

Mdthode selon la revendication 1 , caracterisee en ce que Ton procede successivement a une premiere operation 
de relaxation dans le positionnement des noeuds sur chacune des differentes faces triangulees, le peuplement 
en noeuds internes etant effectue a partir des grilles relax6es, et ensuite une deuxieme operation de relaxation 
dans le positionnement de ces noeuds internes. 

Methode selon la revendication prec£dente, caracterisee en ce que Ton realise la premiere operation de relaxation 
par des phases iteratives de relaxation-projection. 

Methode selon Tune des revendications 1 ou 2 ou 3, caracterisee en ce I'on choisit comme direction de projection, 
une direction perpendiculaire a la grille intermediaire. 

Methode selon I'une des revendications precedentes, caracterisee en ce que Ton utilise un mode d'interpolation 
de type transfini. 

• Methode selon Tune des revendications prec6dentes, caracterisee en ce que Ton rearrange le positionnement des 
noeuds des cellules Elementaires creees de facon a obtenir des maillages reguliers. 

Methode selon I'une des revendications 1 a 5, caracterisee en ce que Ton rearrange le positionnement des noeuds 
des cellules elementaires creees de facon a obtenir un mode determine de progression de la tailte des mailles. 

Methode selon I'une des revendications precedentes, caracterisee en ce que Ton procede a un reassemblage de 
tous les blocs. 
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Fig. 6 
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